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FOREWORD 


This interim re port, "Mathematical Model of a Flexible Space Shuttle 
Vehicle," describes the development of a mathematical model suitable for 
controller design. The model was prepared under Contract NAS8-25708, 
"Controller Design Technology for the Space Shuttle Vehicle, " for the 
National Aeronautics and Space Administration, George C. Marshall Space 
Flight Center. 

Dr. S. W. Winder of the Dynamics and Control Division of the Aero -Astro - 
dynamics Laboratory served as technical monitor for the contract. The 
study was performed in the Research Department of the Systems and Research 
Division of Honeywell Inc. Dr. E. E. Yore served as program manager. 

Dr. C. A. Harvey was the principal investigator. Mr. T. Yam performed 
the necessary computer programming and associated numerical analyses. 
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SECTION I 

INTRODUCTION AND SUMMARY 


This report describes the development of a mathematical model of the lateral 
motion of a flexible space shuttle vehicle during ascent. The model was 
developed to perform control system synthesis using stochastic constrained 
optimization techniques. The goals of the control system synthesis are to 
demonstrate the applicability of the techniques and to discover any problems 
peculiar to the flexible nature of a shuttle vehicle. 

The equations of motion are derived in Section II. A brief description of the 
generation of numerical data is given in Section III. Appendixes contain 
explicit definitions and numerical values of trajectory data and coefficients 
appearing in the equations of motion. 

The "piggyback" configuration of the Shuttle was assumed for this study. 

The North American Rockwell phase B study vehicle identified as the 16 1C- 
B9U served as the specific example represented numerically. Linearized 
equations of motion are derived based on a nominal trajectory for which 
angle of attack was constrained to be zero. Quantitative data for this tra- 
jectory are given in Appendix D. Aerodynamics associated with flexure was 
computed based on slender-body theory and the assumption of no interference 
between the orbiter and booster. The adoption of this assumption was based 
on a private communication with aerodynamicists at General Dynamics, 

Convair, in which they stated that comparison of wind tunnel data using this 
assumption was good. A modal representation is assumed for gust penetra- 
tion. The gust disturbances are modeled as a stochastic process. 

Three further asumptions were made in developing the model. The first 
assumption dealt with ordering flexure modes consistently along the trajectory. 

It was judged that two flexure mode frequencies crossed each other during the 
interval from maximum dynamic pressure to burn-out. Flexure data at 
additional points would be required for verification of this assumption. The 
second assumption concerns inclusion of "local" aerodynamic terms. These 
terms result from assuming that local normal velocity is the cross -flow 
velocity. Fortunately, these terms were found to be numerically insignificant 
so that effort to verify this assumption was not warranted. The final assumption 
was the incorporation of a gimbal logic. This particular logic provides for the 
two gimbals of each of the 12 engines to be driven by two command signals. 

The two commands yield uncoupled roll and yaw moments. 
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Numerical data for the trajectory, rigid body dynamics, and structural flex- 
ure at five mass loadings served as the starting point for generation of the 
numerical model. In addition to the purely numerical tasks of integration 
and interpolation, the numerical model generation required significant 
"bookkeeping" to assure proper identification and consistent dimensions. 
Finally, "debugging" of the numerical data was required. This debugging 
task seems to lack a formal methodology. Particular techniques used in this 
study are briefly described along with the data generation description in 
Section III. 

From a modeling point of view, the number and distribution of mass points 
for which flexure data were available was insufficient. The question of con- 
sistent mode ordering could not be completely resolved and significant sub- 
jective judgment was required for interpolation. From a controller design 
point of view the flexure data were also lacking. The first mass point cor- 
responds to greater than liftoff weight, and the second point corresponds to 
occurrence of maximum dynamic pressure on the nominal trajectory. At 
least one intermediate point would be desirable. However, the model de- 
veloped subject to these limitations is believed to be sufficiently representative 
for controller synthesis study. 
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SECTION II 

EQUATIONS OF MOTION 


Equations are derived for the lateral motion of a flexible space shuttle 
vehicle during ascent. Two independent assumptions made in the derivation 
are ( 1 ) aerodynamic forces on the orbiter and booster are independent of each 
other and may be computed using slender-body theory (References 1 and 2), 
and ( 2 ) side gust penetration may be represented in terms of three "mode 
shapes". The rolling gust is "distributed" over the wing chord with a single 
mode. 

We begin with a discussion of relevant coordinate systems. Then we describe 
our flexure representation. Linearized equations of motion are then derived 
using slender-body theory to represent aerodynamic forces. 


COORDINATE SYSTEMS 

Two body-coordinate systems shown in Figure 1 will be used in this derivation. 
The (xj,yj, z^) coordinate frame (1-frame) is the No. 1 coordinate system 
used in the structural mode computer listing provided. The (xg» y^, ZJ 3 ) co ' 
ordinate frame (B -frame) is a conventionally oriented aircraft body coordinate 
system with origin at the center of mass, which is located at (xj, 0 , z^) in 

the 1 -frame. The vector r(t; xj, y^, z-^) depicts the vector from the origin 
of the earth coordinate system to the point with coordinates (xj, y^, z^) in 

the 1 -frame. As in Reference 3, a pitch ( 0 ), roll (0), yaw (i//) system of 

Euler angles is used to relate the B-frame to the earth frame. Thus, B-frame 

base vectors (L,, 1 -,, k_) are related to earth-frame base vectors (i , i , k ) 

B J B B e J e e 

by: 


l B 


(cdcty + S0S0S1//) (c0C^) ("S0C^ + C0S0S^) 



% 

= 

(-cdsxf/ + SdsQclp) (c0cl//) (S0S»// + C0S0CJ//) 


h 

.V 


(s 0C0) (-S0) (c 0C0) 


k- 


Note that ig = -ij, j B = jj, and k B = -kj. 

The rate of angular rotation, to, of the B-frame relative to the earth frame is: 
W = P 1 B + q j B + r k B 
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Figure 1. Vehicle Geometry and Coordinate Frames 

The rates p, q, and r are related to 0, 0, and i// by: 

p ci/' (c0)(si//) 0 1 p"l 1*0] [“ cl// -si// 0 ] Tp 

q = -si// (c0)(cl//) 0 0 0 = (sl///c0) (cl///c0) 0 q 

, r . .0 “ s0 1_ _ >/^ , . )//. _(tan0)(sl//) (tan0)(ci//) 1_ _r 

The velocity of the center of mass with respect to the earth frame is: 

^(r> = u i B + v j B + w k B 

where r = r(t; 0, 0, 0) - ig - z ^ kg. It will be assumed that u, w, and q 

are not disturbed from their nominal values of u , w , and q. Also the nomi- 
nal values of v, 0 and i// are zero. o o ^o 
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FLEXURE REPRESENTATION 

The flexure is represented by a set of natural vibration modes determined 
from a stick model of the mated vehicles. These modes are defined as lateral 
oscillations of center lines of the orbiter and booster. The reference center 
line of the booster is the xj axis. The orbiter reference center line is the line 
z i = Z 10> yi = 0. Thus, the mode shapes are given as functions Y i (x^, z^) 
defining normalized deflection of the points (x^, 0^ z^) on the center lines to 

the points Cx x , Yj (x 1 , Zj), z 1 ] for discrete values of x 1 and the two values of 
0 and z 1Q for z y Smooth functions Y^Xj, Zj) with respect to x x and partials 
of Yj(xj , z x ) with respect to Xj as needed are generated with spline fits to 
the discrete data. The total flexure is then represented by: 

Y(t; Xj.Zj) =Sti. (t) Y. (x 1 ,z 1 ) 

where r\ 1 (t) denotes the magnitude of the deflection of the i-th mode. A 
typical mode shape is sketched in Figure 2. 



Figure 2.. A Mode Shape 


The normal modes of vibration and their corresponding frequencies are func- 
tions of the mass distribution and structural stiffness of the vehicles. The 
mass distribution changes as fuel is expended and hence is a function of time 
along the trajectory. Thus, in fact, the mode shapes are functions of time 
although we suppress this dependence on t in all of our notation. The mode 
shapes are considered as "slowly varying" with time. This requirement 
brings about the need for special consideration of consistent ordering of the 
modes and uniform normalization of each mode shape all along the trajectory. 
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The flexure data provided for this study consists of modal frequencies and 
normalized displacements at a finite set of vehicle stations computed for the 
161C-B9U using NASTRAN* for five distinct mass distributions. These 
mass distributions correspond to approximate times of liftoff, burnout, and 
intermediate times along the trajectory of 70, 90, and 110 seconds after lift- 
off. The model consisted of 29 grid points, 27 connecting bars, 18 concentra- 
ted masses, and two scalar spring connections. 


Although common grid points and structural properties were used, the modal 
data (frequencies and mode shapes) were computed independently for each 
mass distribution. In each case the NASTRAN output orders modes by 
increasing frequency, and the mode shapes are normalized such that the 
maximum displacement is of unity magnitude. 


The lowest four frequencies are shown in Figure 3 for the five mass distri- 
butions considered. We are faced with the task of constructing four curves 
which represent the frequencies as functions of the mass loadings. The lowest 
and highest frequencies present no apparent problems. However, there is 
some question for the middle two frequencies, especially for M„ ^ M ^ M . 

That is, there is reasonable doubt as to the correct connections of the points 
A and B with C and D. The modal masses corresponding to these two sets of 
frequencies are shown in Figure 4. This appears to be of little help in deter- 
mining the proper connections. Thus, we turn to the corresponding mode 
shapes. These are shown in Figures 5 to 8. For completeness, the mode 
shapes for the lowest and highest frequencies are shown in Figures 9 and 10. 
The mode shape corresponding to point D appears to be much more likely a 
successor to point A than point B. So the connections from A to D and from 
B to C were made. To be certain of the proper connection, data for inter- 
mediate mass loadings would be required. 


Figures 5 and 8 show normalized shapes of what we call the same mode. But 
the point on the vehicle at which the shapes are normalized is the orbiter nose 
for the first four mass loadings and is the nose of the booster at the fifth mass 
loading. The normalization was changed so that this mode shape had a unit 
displacement at the nose of the orbiter for all mass loadings. In this way the 
corresponding r\ (t) may be interpreted directly as displacement of the orbiter 

ln 4 V mts wl \ ich are constant along the trajectory. Similar normalization 
oi the other modes was made where necessary. 


"NASTRAN is a general-purpose digital computer program designed to 
analyze the behavior of elastic structures under a range of loading 
conditions using a finite-element displacement method approach. 
Reference 6. 
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Mass (millions of pounds) 

Figure 3. Bending Frequencies for the 161C-B9U 
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Normalized deflection 


® Orbiter deflection 
a Booster deflection 

Mode shape at mass = 5.2 x 10 6 lb. 

Mode shape at mass = 3.9 x 10 6 lb. 

Mode shape at mass = 3.6 x 10 6 lb. 
Mode shape at mass = 3.2 x 10 6 lb. 


_ . - 
oo* ^ 3 *' 


/ M / 


/ // 
s // 

7 / 


/ i 

S / 


Distance from nose of orbiter (in.) 


Figure 5. Second Mode Shapes at First Four Mass Points 



Normalized deflection 


° Orbiter deflection 

* Booster deflection 

Mode shape at mass = 5.2 x 10^ lb. 

Mode shape at mass = 3.9 x 10^ lb. 

Mode shape at mas$= 3.6 x 10^ lb. 

Mode shape at mass = 3.2 x 10^ lb. 



Figure 6. Third Mode Shapes at First Four Mass Points 
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Normalized deflection 


5 


o Orbiter deflection 
a Booster deflection 


a 
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Figure 10. Fourth Mode at Five Mass Points 


LINEARIZED EQUATIONS OF MOTION 
The linearized kinematic equations are: 


0 = p - Q q ^, ^ = r + Q q 0 , y e = v + u Q tf/ - w Q 0 


The side-force equation is: 


-f * 3 B > = -f" ^ + u o r " w o p) 

N M 

= W[0cos 6 sin0 ] + ]• *[S F„ +S F 5 + jjF.dS] (1) 

° ° B 3 = 1 ] 3 = 1 3 S 
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Where Ft- = thrust of the j-th engine, Fg^ = j-th control force. Fa = local 
aerodynamic pressure, and S denotes the surface of the vehicle. 

The j-th engine is located at (x-^g, yE B zjr;. ) the 1-frame. The thrust of 

J J 

each engine is assumed to have magnitude, T/N and to be directed forward 
parallel to the "deflected" booster reference line. Variations due to gimbal 
deflections are contained in the control forces. Thus, 


T. 

3 


T 

N 


Ci 


B 


3B SY i' (x lS' 0), li ] 


The j-th control force is assumed to act at the point (xg., yg , zg ). With the 

thrust and control forces acting at these specified points, the mofkent equations 
about the x B and z B axes are: 


r 


— — 


r n 


- - 


— — 

i 

-I 


p 


i 

I - I 


p 


0 

XX 

X z 




xz 

yy zz 


Jr 


-I 

I 


r 

= Q 

I - 1 

-i 


r 

+ 

M. R. 

xz 

zz 




XX yy 

xz 





L_ 


— — 


L J 


— « 


— 


N 

+ S 

3=1 



j B <z E.' z l> +k B SY i (x l6 ,0)r,. 
3 

J B (x l " x l6 ) ~ i B SY i (x lg’°> T l i 


M 

+ S 

3=1 



j B (z 6. ■ z 2 ) +k B EY.(xg ,0 )ti. 
3 3 

j B ( Xi - Xg ) - i B EY.(xg ,0)ri. 

3 3 


+ 



(z i ■ z i ) j B + yi k B 

(x r x i)j B - yi i B 


F A d S 


( 2 ) 


(3) 
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where M. R. denotes the inertial rate and jet damping term of Reference 3. 
It is assumed that the x-^ -plane is a plane of symmetry. The equations 
for the mode deflections are* 


m. ffj.+o. 

1 '1 




+ w, 2 T 


) = 3 


B 


N 

•E 

3=1 


Y i< x l6’ 0) 


+ 


'B 


M 

£ F 
3 = 1 


6 Y i (x 5.' 0) 


+ ^B F A dS 

The thrust contribution to Equation (3+i) will be neglected since the thrust 
is assumed always to be normal to each of the mode shapes. The thrust 
contribution to Equations (1), (2), and (3) are: 


-t£ y.' (x i6 ,o)ti. 
i 


T_ 

N 


•S [Y i '(x l5 ,0)r|. S (z E - 

3 3 

- TS Ux 1 - x) Y. # (x l6 ,0) +Y.(x i6 ,0)] R. 


The control force contributions to the equations are taken directly from 
Reference 3 for Equations (1), (2), and (3). They are: 


qSC 6a + T 
J 6a 


qSBC. 6a + T 
^ 6a 


qSbC 6a + T 

n 6 a 


L \ 


^l«p + 1-^- 


'6p 


6r 


N 


«E_ 


j 6p + 

( L 6r ' 

6r 

T 


i ; 

- 

6p + 

1 N 6r ] 

6r 

1 T 

i 1 

1 i 

- 


*A structural damping term (0.2(0/^) was arbitrarily hypothesized. 


(3+i) 

(la) 
(2 a) 
(3 a) 

(lb) 
(2 b) 
(3b) 
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In these equations 6a denotes aileron deflection. The quantities 6p and 6r 
denote roll and yaw control gimbal deflections defined as follows. Each engine 
is assumed to be gimbaled such that 6yj and 6zj represent positive rotations 

of the j-th engine about axes parallel to the y B and z B axes which are centered 
at (xjg , y E , z E ). The gimbal deflections of the j-th engine are defined as: 


6y . = c 6p 
yp-j 


6z. 

3 



6p + c 


zr. 

3 


6r 


1 x 

— sgn(y E .) where sgn (0) = 0 and sgn (x) = — 

3 M 


if x t 


Engines of the B9U are arranged as shown in Figure 11 in which the positive 

xi axis is directed out of the paper. The magnitudes of c n . and c are 

Pj zr- 

constrained not to exceed unity. The values of the coefficients are determined 


so that maximum rolling moment per unit 6p is to be produced subject to the 
constraint that 6p produces no yawing moment about the velocity vector, and, 
similarly, maximum yawing moment per unit 6r is to be produced subject to 
the constraint that 6r produces no rolling moment about the velocity vector. 



Figure 11. Engine /Vehicle Geometry 
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The control force contribution to the mode deflection 


T Y.(x. 6 , 0) 


6p 


5p + 


6r 


equation (3+i) is: 
C(3+i)b] 


The remaining contributions required are those from the aerodynamic force 
F A . It is assumed that the aerodynamic force may be partitioned into two 
parts, one acting on the orbiter and one acting on the booster. This may be 
expressed mathematically by setting 


F A (t; x 1 , Zj) = F Q (t; Xj) 6(z 1 


— z 1 q) + F b^' x l^ 6 ^ z l^ 


with the definition thatjjg F(t; x, z) 6(z - z Q )dS = F(t; x, z^) dx. 

Slender-body theory yields the following expressions for the local aerodynamic 
side forces per unit length: 

nq f~ . o dv o< x l> t} 

P o <t - x l ) • % * 2V f “oW <*1> v o<*r‘> + D 0 (Xj) 8x * ■ 


. 1 - 2 . . 

+ ^ D (x, ) 


V ~o 


at 


9 v b( x l't) 


F b <t ' x l > ' % = ^j 2D b (x l )D b (x l )v b (x l' t) +D b <x l> S Xl 


+ i-n 2 
V D b (x^) 


3v o (xi,t)| 

at 


where D Q (x^) is the cross-sectional dimension of the projection of the 
orbiter onto the x^z^ -plane and v 0 (xj,t) denotes the j^-component of the air 
velocity with respect to the orbiter. and v^(x^,t) denote similar 

quantities related to the booster. Thus, 

v o (x r t) = -r (t; x r o, z 1Q ) • jg + v w (t, x^ 
v b (x 1 ,t) = -r (t; x 1 , 0, 0) • J B + v w (t, x 1 ) 
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The local cross-wind, v w (t, x^), is assumed to be independent of Zj and 
represented by penetration states: 

3 

v w (t ’ x i> = E 

i=l 

where 

f-^x) = 1 - exp (-2.3 x/L) 

f 2 (x) ~ 1 _ C 1 - f^( x)] cos (2tt x/L) 

fg(x) = f 2 (x) - [1 - f^x)] sin (2tt x /L) 


and the xjt) are states driven by the cross-wind at the nose of the vehicle. 
That is: 


x„ 



v 

w 
V , 




V_ 

L 


3.99 X2 - 6.29 x 0 + 2.30 -rr 


w 


V 


x = ^ 
x 3 L 


12.58 x 2 - 8. 59 x 3 - 3. 99 ^ 


with v w - v w (t’ x iisr)- The choice of the c^s will be described later. The 
wind v w (t< x in^ is assumed to be (Reference 5). 

v w (t ’ X 1N } = * + CT w u 

where 


r • "i 






p 

, to 


0 h 


to 


1.378 * 10 -2 ' 


= 

-c c h -c.h 

L 5 4 J 


X 


-0. 965 • 10" 6 


with c 4 = 1.9 ' 10 -4 , c 5 = 1.44 • 10' 8 , 

CT W = standard deviation of the wind 

h = altitude 

n. = unity white noise 


18 



The position vector, r(t; x^, 0, z^), is given by 

rftjXj, 0, Zj) = Htjxy 0, z j) + (x 1 - Xj) ig + (z - z-^kg 

+ ^B ?' VVV n i 

where r(t; Xj, 0, z^) = r(t) and r (t) = u Q ig + v jg + w Q kg • Thus 
r(t; Xj.O.Zj) = u 0 i B +vj B +w 0 k B +S YjSj, j B +<*!’ 

+ (Zi- Zl )kB +j B f Y i <x l- Z l )r lr 

Using the relation (linearization being assumed) 


_d 

dt 


a B 


0 r -Q„ 


-r 0 p 

Q -p 0 

o 



l B 


^B 


k„ 


B 


leads to 


rijtjXj.O, Zj^) = [u Q + (z 1 -z 1 )Q Q ]ig 


+ [v + (x^-x^r - (z^-z^p +S Y.fxj, z 1 ) T l^j B 
+ [w Q - (x 1 -x 1 )Q o ] kg 

Thus r(t; x y 0, zj * jg * v+ -(Xj. - x^r -(^-z^p +S Y.(x r 


and 


^ [r (t; x v 0, Zj) • jg] = v + (xj-x^r - (z^z^p +S Yjtej, z^r] 


The results may be combined to write: 



2 V 

F 

nq B 


( t, x) 


j B = 2D b (x)D b , (x) [-v - (x 1 ~x)r + ZjP -E Y.(x, 0)^ 

i 

3 3 

+ E c.f.lxlxj] +D b (x)[r - S Y. (x, 0)f|^ + S c,f/(x)x.] 
i=l i i=l 1 


+ Y"D b (x) [-v -(x-x)r + z.,p -Ey.(x, 0)r). +E c . f.(x)x.] 
* ^ * i * i li i 

A similar expression can be written for the orbiter. These forces may be 
entered into the integrands in Equations (1), (2), (3), and (3+i). 

In this derivation, v b (x, t) and v Q (x, t) do not contain the terms: 

u o? Y iV o)T1 i and u o? Y i' (x l' z 10 )r V 


These terms appear in Reference 2 and have been included in previous studies. 
(References 4, 5). This type of term would arise if the local normal velocity 
is used in the definition of the cross-flow. That is, defining for the booster 

v b (x r t) = -«t; XjO.O) • S Y.'lxj.o) nj] +v w 


and then consider ( z.Qq) negligible relative tou Q . Suchu Q terms were 
arbitrarily added for the booster and orbiter in the following expansion of 
equations (1) - (3+i). These terms are underlined. Fortunately, the numerical 
magnitudes of these terms did not significantly modify the system in the sense 
that eigenvalues computed with the terms differed insignificantly from eigen- 
values computed without these terms. Resolution of the correct treatment of 
these terms was beyond the scope of this study. 


The following form of Equations (1) - (3+i) were obtained by making the 
appropriate substitutions and using integrals defined in Appendix A. 


Equation (1) is: 
v 


W + H 
g V 1 ! 


+ | (V + V + . 2 % \ 


6+1 T 1 i 


W 


= P(w X 


of ' 2HI 11> +r[ - u oI +H(I r 2I 12 )] - 2HI 13 v+2 WljHncoso -l)I 13+i 


i=l 


■ 2I 17 + i ] - TY i ,(x l6'°> J + i 2x i Hc i< I 29 + i +2I 32 + i ) 


+ W(0 cos 0 Q + xjs sin 0 Q ) + q S 5a + T 

where H = -^-2. and four flexure modes are included. 
’ 2 V 


17%) 

f Y « r l 


6p + 

\ T / 

i T i 

6 r 


20 



Equation (2) is 
H 


<‘xx + t ! 36>P - «xz - § V - 2 Ciiil37+i + El "* l40+l1 


= P(Q 0 I XZ + 2 HI 45 ) +rCQ o < Vy"W " Ha 2' 2I 46 )] t2HI ll T 

+ ^ {^HC cos g 0 “ 1 ) I 46+j ‘ 2 I 50+i^ + T| i H u o^ 54 +i + 2 I 58 +i^ 


+S X, H c,(Ljo , t + 21^44) +q Sb C. 6a + T 


i=l 
L 


+ 


1 - i' 62 +i 65 +i 

N 

£ 

N ]=1 “j 


6a 


J 6pl 


6p 


L\ 


6r 


Equation ( 3 ) is: 


ar - l s , < Z E, -V S =1 Y i'< x l6-°>\ 


3 . 3 


(I zz + V "69 


f -Ov, ‘I + 7 C V 'E ^ l69+i + ?=l’ 11 l72+ ‘ ] 


xz V 37 F V 3 i=1 

PKU'xx'W ‘ 2H X 46 ] + r[ - Q o'xz + H<1 3 ' 2 ! 77 )] ' 2H *12 v 

4 
£ 
i=l 

- Wij-XjjJYj'fXjj.O) - T Y.( Xl g, 0)] + S Xj He, < 93+i + 2I 96+i > 

N._\ 


*0 xx yy‘ 

4 

+ £ ^.H[ (cos a Q - l)I 77+i - 2Ig 1+i ] + ^' 1 T 'i^ H - u o^ 85 +i + ^ 89 +i^ 


Ic I l>lf 

+ q Sb C 6a + (M.R.) r+ T -£&) 6 P + r 
4 n 6a T / IT 



6r 
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Equation (3+j) for j = 1, 2, 3, 4 is 
4 


E 6 ij M j + y + y ^40+jP + + I 6+j' 


-S c.x I (i, j)] 
1 = 1 1 1 * 


=2p H I 50+j + r H(I 64j - 2I 81+j ) - 2Hv 1^ - M. [0. 02^ . 
+ W ?T1.] +S .H[(cos ot - 1) I 3 (i, j) - 2I 4 (i,j)] 


3 3 


i=l 


+ T!. Hu Q [l 5 (i,j) + 2 Ig (i, j )]} 


3 

+ s 


P x i U c. [I 7 (i,j) +2 Ig(i, j)] + T Yj( Xl g, 0) 


L\ 


6p 


6p + 


Some of the integrals in the above equations represent familiar terms. 
For example, 

I, 



Sb 

11 = 

" l /3 


Sb _ 

12 = 

'~S 


s ^ 

13 = 

--r C 

n yp 


with the center of mass being the moment reference points for C, and C . 

n j3 

These integrals will be replaced by the wind tunnel data. Noting that for a 
constant cross-wind input, say v (t) = v w , the steadystate value of each of 
the xj's is v w , and observing that in this case the wind input to Equations (1), 
(2), and (3) are: 
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H i?i X i C i (I 29+i +2I 32+i ) V i?l c i (I 29+i + 2 I 32+i ) 


H S x i c i ( 62+i + 2I 6 5+i ) 
i=l 


H 2 x.c. (I OQ . . + 2I ne .. ) 
i = j 1 1 93+x 96+i' 


w 

V 


w 

V 


[ V 2 , c i a 62 + i + 2 W 
1 = 1 

*^J S c i (I 93+i + 2I 96+i ) 

'Vi 


we may satisfy correct steady state side-force, yawing, and rolling moment 
inputs by choosing the c^’s as solutions of: 


(I 29 + i +2 W c i SC 


y /3 


- 3 


i3S (I 62+i + 2I 65+i )c i =qSb % 

2 i = l P 

— r (Iqq+s + 2I oc ..)c. = q Sb C 
„ Zj 93+1 96+1 l n no 

1 i=l P 

The c^'s will be so chosen with Cy , C^, and C n ^ representing wind tunnel 
data with the center of mass being the moment reference point. 

The rolling gust is represented by (Reference 3) 

P g = d 24,24 P g + f 24, 3 n 2 

where n 2 = unity white noise. The rolling gust is "distributied" by X 4 over 

the wing chord with steady-state values of x. and p being identical. Thus, 

4 g 

X 4 = d 2 3 , 2 3 X 4 + d 23, 24 P g 

with dgg 23 = -< ^23 24’ state > x 4 * drives the vehicle in the same manner 

as geometric roll rate p. Hence we may substitute p + x for p in the right 
hand sides of Equations (1), (2), (3), (3+j). 
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SECTION III 

GENERATION OF NUMERICAL COEFFICIENTS 


The equations of motion derived above may be viewed in vector form as 
Cx = Dx + Eu + Fd 

where C, D, E, and F are matrices of appropriate dimensions and 

^ ” CP* Y » 'Hj > Tlj , Tig » Tig j "Hg > T|g , r^, Tjj, 5p, fir, 6a, X ^ ,Xg , Xg 

x,u,x 4 ,p g ] ' 

u - [u r u 2 , u 3 ]' 
d = [ v, n , n 2 ] 

These equations of motion describe perturbation about a nominal trajectory 
which is assumed to be defined with time as the independent variable. Thus, 
for a given trajectory, the matrices C, D, E, and F have continuously time- 
varying elements. For controller synthesis, values of these elements were 
computed at endpoints of five-second intervals along the nominal trajectory. 
These data were then used to compute the coefficient matrices at the same 
time points for the equations of motion in the form 

x = Ax + BjU + BgV + Bgn 

with n - [n , n„] . These coefficients, along with those for the response 
equation 


r = HjX + Du 

are listed in Appendix C. The manner in which these coefficients were com 
puted will now be summarized. 


The original data consisted of a nominal trajectory, tabulated moments of 
inertia and center of mass location at four mass points, tabulated rigid body 
aero data such as Clp at 10 values of Mach number, and tabulated mode de- 
flections at discrete points of the booster and orbiter for five mass points. 

of Parties and rigid body aero data at intermediate points 

corresponding to the five-second time points on the trajectory were ob- 
tained by linear interpolation. 
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Integrals over the vehicle as listed in Appendix A were computed using cubic 
spline-fits (Reference 7) to the discrete representation of the mode shapes and 
cross-sectional diameters. The integrals involving mode shapes or their 
derivatives were evaluated at the five given mass points. The resulting inte- 
grals were then spline-fit to give values at intermediate mass points corre- 
sponding to the "five- second" time points. Similar spline- fits were used for 
all other coefficients involving modal data. The five mass points provided 
were a minimal number to base the fits on with expectations of confidence in 
the resulting functions. The fits were quite sensitive to assumed end-point 
slope conditions. Subjective judgment played a significant role in choosing 
these conditions. Data for two more points, one near each of the end-points, 
would have been quite useful in determining these boundary conditions. 


After evaluating the necessary integrals, the coefficient matrices C, D, E, 
and F and, in turn. A, Bi, B 2 , and B 3 were computed for the 43 time points: 
0, 5, 10, . . . , 210 seconds. At each of these time points, a sampled-data 
version of the system was computed. The sample-data version is of the form 


k +1 


= A 


x k + 


u k + 


V k + 


B, 


“k 


where 


A = exp [A(At)] 

B. = A _1 (A - I)B. 


In the right-hand side of 


At [J + 


(A At) 


+ 


the last equation, A 

(A At ) 2 . (AAt) n 

3! • ‘ ‘ + ?n+TT! 


(A - 
+ . . 


I) denotes the series 

] 


The A, B^ matrices were computed using the series representations and 
truncating the series after including a sufficient number of terms to guaran- 
tee convergence to at least three-digit accuracy. Controller synthesis and 
analysis was performed using these coefficient matrices with linear inter- 
polation for intermediate points. 

Verification of the numerical model is an extremely difficult task. Three 
specific techniques used in this study were examination of time-varying coef- 
ficients for unusual algebraic signs, sign changes, and magnitude variations; 
eigenvalue comparisons; and response examinations. The first method was 
useful in finding errors resulting from some keypunch errors, certain alge- 
braic errors in equation derivation or similar errors in programming, and 
unreasonable effects of assumed slope conditions in spline -fits. The second 
method consisted of comparing eigenvalues of the A matrix with known rigid 
body characteristics of similar vehicles and with the uncoupled flexure char- 
acteristics. This method proved useful in finding algebraic and programming 
errors. The third method consisted of determining whether computed responses 
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could reasonably be expected from the inputs used. This method led to dis- 
covery of keypunch, algebraic, programming, and dimensional errors. Un- 
fortunately, human errors do occur, and at this point model verification is a 
tedious and somewhat unorganized task. Development of standard, objective, 
and direct techniques for model verification would be of great benefit. 
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APPENDIX A 
TABLE OF INTEGRALS 
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APPENDIX A 
TABLE OF INTEGRALS 


The lower and upper limits of integration are x 1N and x 1T for each of the 
following integrals. 1 

h = J [d 0 2(x) + D b 2(x) ^ dx 

1 2 = J [(z r z 10 )D o 2(x) + z 1 D b 2 (x)]dx 

13 = J (Xj-x) [D q 2 (x) 4- D^ 2 (x)] dx 


I 3+i 

= -J f.(x)[D o (x) + D b ^(x)] dx 

i = 1,2,3 

I 6+i 

= S [Y.(x, Zl0 )D o 2 (x) + Y.(x,0)D b 2 (x)]dx 

i ■ 1,2, 3, 4 

'll 

= J [ <V Z 10 )D o (x) D o ,(X) + z i D b (x) D b /(x)] dx 


*12 

= J* (x 1 -x) [D q (x) D q '(x) + D b (x) D b ' (x)] dx 


'l3 

= J ED o (x) D q '(x) + D b (x) D b '(x)] dx 


'l3+i 

= J ^D o 2 (x) Y. ' (x, z 1Q ) + D b 2 (x) Y. ' (x, O)] dx 

i = 1,2, 3, 4 

*17+i 

= J[D q (x) D q '(x) Y.(x, z 1Q ) + D b (x) D b '(x) Y.(x, O)] dx 

i = 1,2, 3, 4 

J 21+i 

= J [D q 2 (x) Y. "(x. z 1Q ) + D b 2 (x) Y. ' (x, O)] dx 

i = 1,2, 3, 4 

J 25+i 

= J* [D 0 ( x )D o /( x)Yi '(x, z 1Q ) + D b (x) D b '(x)Y. '(x,0)] dx 

i = 1,2, 3, 4 

*29+i 

= J V (x) [D q 2 (x) + D b 2 (x)] dx 

i = 1,2,3 

: 32+i 

= J f.(x) [D q (x) D q / (x) + D b (x) D b 7 (x)] dx 

i = 1,2,3 

'36 

= I t(z r z 10 ) 2 D q 2 (x) + z x 2 D b 2 (x)] dx 


*37 

= J <x r x) CD o 2 (x) - z x D b 2 (x)] dx 


J 37+i 

= J* fj(x) [D q 2 (x) (z 10 -i x ) - Zj D b 2 (x)] dt 

i = 1,2,3 

' 40 +i 

= J* [Y.(x, z 10 )D q 2 (x) (z^-Zj) - Y.(x, O) i 1 D b 2 (x)] dx 

i = 1.2, 3, 4 
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J [ <V Z 10 )2 D o (x) D o' (x) + 5 1 D b (x) D b /(x)] dx 
J* (x r x) [(z r z 10 )D o (x) D o / (x) + z x D b (x) D b '(x)] dx 


45 


46 


l 46+i * J[<- 10 - S l )Y i' (x - Z 10> D o 2<x) - dx *■ l - 2 ' 3 - 4 

J[z -5 1 )Y i (x,z 1 0>D o (x )D o'W - ijYjfx.O) D b <x)D b '(x)] dx 

i ■ 1,2, 3, 4 

J[(x 10 -i)Y 1 " ( x,z 10 )D o 2 (x,-i 1 Y l ' (x.O)D b 2 (x)] dx i- 1,2, 3. 4 
J[(z 10 -i)Y i '(x,z 10 )D o (x)D o '(x) - z 1 Y 1 '(x,0)D b (x)D b '(x)]dx ^ j 4 

J fj'(x) K* 10 -i l )D o 2 (x) - ij D b 2 (x)] dx i = 1,2.3 

*65+1 * J" f j (x>[(z l0" z l )D o (x)D o (x) ‘ i l D bW D b' Wldt 1 = l ’ 2 ' 3 


59+i 


l 54+i 


58+i 


62+i 


69 


6 9+i 


l 72+i 


J (x 1 - x) 2 [D q 2 (x) + D b 2 (x)] dx 

J L(x) (x^x) [D q 2 (x) + D b 2 (x)] dx 

f(x f x) CY.(x,z 10 ) D o 2 (x) + Y.(x, O) D b 2 (x)3 dx 


i = 1, 2, 3 
i = 1,2, 3, 4 


I 77 = J(x r x) 2 [D q (x) D q '(x) + D b (x) D b >)] dx 

I 77 .fi = J (Xj-x) [Y. '(x, z 10 ) D q 2 (x) + Y. # (x,0) (x)] dx i = 1, 2,3,4 

= J(x r x) CY.(x,z 10 )D o (x) D o / (x) + Y.(x,0)D b (x)D b / (x)]dx = 

= J (x r x) [Y. // (x,z l0 )D o 2 (x) + Y.'(x,0) D b 2 (x)]dx i= 1,2,3,4 

= J* (x x -x) [Y i '(x,z 10 )D o (x)D o / (x) + Y i / (x,O)D b (x)D b / (x)]dx i = 1,2, 3, 4 

io 3+i • / f i '< x) (x r x) [d o 2(x) + D b" (x)] dx 


81+i 


85+i 


89+i 


I = Pf.(x) (x -x) CD (x) D '(x) + D b (x) D'(x)] dx 

QRj_i 1 1 O O U V 


i = 1,2,3 
i = 1,2,3 


96+i 
I j(i.3) 


- J , ^ i (x,z 10 )Y j (x,z 10 )D o 2 (x) + Y.(x,0)Y ;i <x,0)D b 2 (x)] dx _ ^ ^ 


yt,j)-Jf 1 WCTf j <x.* 10 )D o ; 2 (x) + Y j (x,° > D b (x)] dx 


i- 1,2,3 
j = 1, 2, 3, 4 
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I 3 (iJ)= J* [Y. 7 (x, z 1Q )Y (x, z 1Q ) D q 2 (x) + Y. '(x.O) Y (x, O) D 2 (x)] dx 

i» 3 = 1.2, 3,4 

I^(i» j) = J [Y.(x, ZjqJY^x, z 10 )D o (x)D o / (x)+ Y.(x,0) Y..(x, O) D b (x)D b '(x)] dx 

i,j = 1,2, 3, 4 

I 5 (iJ) = J‘[Y. // (x, z 1Q ) Yj(x, z 10 )D q 2 (x) + Y. "(x,0)Yj(x,0)D b (x)]dx i,j = 1,2, 3, 4 

I 6 (i. 3) =J CY. '(x, z 10 )Y.(x, z 10 )D q (x) D o ' (x) + Y. '(x, O) Yj x, 0)D b (x)D b '(x)] dx 

i, 3 = 1, 2, 3, 4 

I 7 (i,j) =Jf. , (x)[Y ; .(x, z 1Q ) D q 2 (x) + Y..(x, 0) D b 2 (x)] dx i = 1,2,3; 

j = 1, 2, 3, 4 

I 8 (i 'j ) = I f i (x) [Y j (x ' z iO )D o (x) D 0 /( x) + Y.(x.O)D b (x) D fe / (x)] dx 

i = 1,2,3; 
j = 1,2. 3, 4 
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APPENDIX B 

EXPLICIT REPRESENTATION OF COEFFICIENT MATRICES 
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APPENDIX B 

EXPLICIT REPRESENTATION OF COEFFICIENT MATRICES 


The equations of motion may be written in the form; 
Cx = Dx + Eu + Fd 


where C, D, E and F are matrices of appropriate dimensions 

x = Cp, r, v,0, \Jy, n^rig, ri 2 , ri 3 , ri 4 , r\^, 5p, 5r, 5a, Xj, x 2 , x g , x, u, x 4 ,p ] ' 

o 

U = [Uj, Ug, Ug] ' 
d = [v, n p n 2 ] ' 


The matrices C, D, E, and F are actually quite sparse. In fact, the only 

nonzero elements of E and F are, e 15 _ i ei5i2 , e, 7 f, 8 ,. f 19 ,, f 20 , 


! » 2». I 24, 3‘ 


T’U X ■ * „ *■*,*’ VU’ *18, J.» l 2U,l« 

The matrix C has 17 rows and 14 columns of an identity 

t. n ri n A 44 . . ** 

2 . 


^21 2‘ f 22 „ 

matrix, that is7"cij = 6ij for i ^ 1, 2, 3, 7, 9, 11, 13, and Cii = 0 for i = 1 
3, 7, 9, 11, 13 and j = 4, 5, 6, 8, 10 12, 14, 15, 16, 17, 21, 22, 23, 24. 
In the D matrix there are four rows with a single nonzero element, namely 

d . . — J 5 . , . - . C* J _ o i n i a nir • , , . _ * J 

1 1 - 


ij = ®i(j+l) ^ or f _ 10, 12, 14. Major blocks of zeros in D are: 


djj = 0, i ? 3, 4, 5, 6; j = 4, 5, 6 


d ti - °- 


i < 18, j = 21, 22, 24 


dy = 0, i = 4, 5, 6, j > 5 
d ij = 0, i > 14, j < 15 

The lower right-hand (10 x 10) submatrix of D is nearly diagonal containing 
only 17 nonzero elements. The nonzero elements of matrices C D E and 
F are defined explicitly as follows: 


C 44 = C 55 ~ C 66 = c 88 = C 10, 10 = C 12, 12 = 1 


c.. 

li 


= 1, is 14 


'll 


= 1 + S 

xx V 


I 


36 


32 



= T + ii T 
C 12 ■ xz V X 37 


- ii t 

c 13 ’ V 2 


H 


1,5 + 2j V *40 + j, J = lr 2 ‘ 3 ’ 4 


c l, 17 + j V c j I 37+]‘ * 3 - 2, 3 


c 2 1 ~ C 12 


- T + ii T 

c 22 " zz V 69 


H j 

c 23 “ V 3 


c 2, 5 + 2j = V X 72 + j, 3 = 1* 2, 3, 4 


H 


, 17 + j = " V c j *69 + j, 3 = 2, 3 


c 3 1 = ~ C 1 3 


c 32 = c 23 


_ _W , _H_ j 
c 33 g V h 

c 3, 5 + 2j = “ *6 + j* $ ~ ll 2> 3> 4 


H 


c 3, 17 + j " V Cj 


j h + r i- 1 - 2 - 3 


H 


+ 2i, 1 ~ V *40 + i’ 1 = lf 2 ' 3 ’ 4 


c 5 + 2i, 2 V I 72 + i’ 1 2 ’ 3 ’ 4 
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i = 1, 2, 3, 4 


c 

c 

c 

d 

d 

d 

d 

d 


- H T 

5 + 2i, 3 V x 6 + i' 


H 


5 + 2i, 5 + 2j = 6 ij M i + V 

5+2 i, 1 7+j = " ~ C j V 3 ' 1 *’ 1 = 1,2,3j 


11 = Q o 1 xz +2HI 45 

12 " Vyy 


= JSb 
13 " V 


' £ 


1, 5+2j ^ COS Q 'o" 1 ^ I 46+j " 2I 50+j^ * 


l,6+2j 


< = Hu (If. . . + 2I co .) -4r Y/ 


o 54+j 58+j N ] 


1,15 



4,17 = ^ Sb 


6a 


d 


1,1 7+j Hc i^62+j + 21 65+j^ 3 = 1,2,3 


d l , 23 = 2HI 45 


d 91 = Q (I -I ) -2H L e 
21 xx yy 46 


d 22 = -Qo I x Z +H(I 3 - 2I 77 ) + M * R ' 


j = 1, 2, 3, 4 
4; j = 1,2, 3 


j = 1, 2, 3, 4 
N 

(x 0) 2 (z -z ),j = l,2,3,4 
0 k=l ^k 1 
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d 23 = C n, 


d 2,5+2j * H[(cos % -l)X 77+r 2I 81+j ], j-1.2,3.4 

d 2, 6-H2J * Hd o«85 + j + 2 W - T(i; r X 16 >Y i' (x 16- °> ‘ T °>' 

j - 1,2, 3, 4 


_6£ 


* 2,15 


d 2, 16 T \ T 


J 2, 17 


= qSbC 


d 2, l*7+j = Hc 3 a 93+j + 2I 96+j>' j = 1 ‘ S ' 3 


d 2, 23 " ' 2HI 46 


, W qSb n 

d 31 = w o I" ' V C l 


= -u 21+ ip c + HI. 

O , g V n„ 1 


i - £§. c 
33 V ^y. 


d D/l = W cos 6 
34 o 


d 35 = w sin e o 
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d 3,5+2j H[(e° s Q ' 0 -1 ) I i3 + j- 2I i? + j] - j 1,2, 3, 4 


d 3, 6+2 j Hu o^21+j + 2I 25+j^ " T Y j / ^ x 16, j - 1- 2, 3, 4 


d 3, 15 ~ T 


d 3, 1 6 T 


d 3, 17 = q S C y 



6a 


d 3, 1 7+j H c j^29+j + 2I 32+j^ j " 2 > 3 

^ _ qSb 


‘3,23 


d 41 = 1 


d. c = -Q 
45 


d 52 = 1 


d,, = Q 
54 o 


d 63 = 1 


d„ . = -w 

64 o 


d„ - = u 
65 o 


V C l 


P 


5+2 i, 1 


2H I 50+i’ 1 = 1 > 2, 3, 4 
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d 5+2i,2 
d 5+2i, 3 
d 5+2i, 5+2 j 

d 5+2i, 6+2j 

d 5+2i, 15 

d 5+2i, 16 
d 5+2i, 17+j 
d 5+2i, 23 
d 15, 15 
d l 6, 16 
d l 7, 17 
d l 8, 18 
d l 8, 22 

d 19, 19 


H ( I 6 + i - 2 W ' i = 1 ’ 2 ’ 3 ' 4 

-2H I 1?+i , i = 1, 2, 3,4 

= H[ (cos » -l)Io(j,i) - 21. (j, i) ] - 0.026..U.M. 

O 0 ^ lj 1 1 

i,j = 1,2,3, 4 

= Hu Q [l 5 (j, i) + 2I 6 (j,i)] -5 ij u i 2 M i i, j = 1, 2, 3, 4 

T Y^(x^ g, 0) | j , i = 1 , 2, 3, 4 

' TY.(x lfl , 0) - i -1,2,3, 4 

= H c.[l 7 (j,i) + 2Ig(j, i)] , i = 1,2, 3,4; j = 1,2,3 

J 

d 5+2i, 1 
-31.6 

-31.6 

-10 

-2.3 V/L 

-<°v /V)d 18, 18 
w 

3.99 V/L 
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-6.29 V/L 


d 19, 20 
d 8 , 7 = 
d 19,22 

d 2 0, 19 
d 2 0 , 2 0 
d 2 0 , 22 

d 22, 21 

d 21, 22 
d 21, 21 
d 23, 23 
d 23, 24 
d 24, 24 
e 15, 1 = 
e l 6 , 2 = 
e 17,3 


d 10, 9 “ d 12 , 11 " d 14, 13 


” (a v / V ^ d ig 19 + d l 9 20 ^ 
w 


= 12. 58 V/L 


= - 8 . 59 V/L 


- (<T v ^ V ^ d 20, 19 + d 20, 20 J 
w 


= h = altitude rate 


= -h [ (0. 95 ) 2 + (0. 735)2] 10 
= - 1 . 9h ( 1 0 ~ 4 ) 

= -2.3 V/(r.c.) 

= -d 23, 23 
= ... -tt V/ 4b 


“ d 15, 15 
-d l 6 , 16 
“ d 17, 17 


1 
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f 19,l = " (d 19, 19 + d 19,20 )/V 

f 2 0, 1 = ”^ d 2 0, 19 + d 2 0, 2 0)/V 

f 21,2 = f 2 2,2 t[(0 - 95)2 +( °- 735>2]1/2 ' 1,9} 10 " 4 
* 22,2 “ 10 ' 2 V 


l 24, 3 


na p 

"5E 



The inertial rate and jet damping term appearing in d 22 is: 


M.R. 


(x 


16-1 + 15)2 

32.17 


d 

dt 



I . 

zz 
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APPENDIX C 

NUMERICAL VALUES OF COEFFICIENT MATRICES 


40 



APPENDIX C 

NUMERICAL VALUES OF COEFFICIENT MATRICES 


In this appendix, the coefficients associated with the equations of motion are 
displayed at five-second intervals, from 0 to 210 seconds. In Appendix B 
the equations of motion were written in the form 

Cx = Dx + Eu + Fd 

Multiplying the equation by C" 1 , and separating d into two components, the 
following standard state vector representation is obtained: 

x = Ax + BjU + BgV + B g n (Cl) 

C _1 D * A 
C _1 E = 

C -1 F = [B 2 B 3 ] 


n = [n^ng] 

The response vector r is defined to be 

r = [0,0, y, y, 5p, 6p, 6r, 6r, 6a, 5a, q/3,q/3, a y , a y ] 
where 

q/3 = [v + au + v ] 

4 

a = v + (x x - 11 0)r - (S 1 - 40)p + £ Y.(110, z 1Q )^. + U Q r - w q p 

y i=i 

-g(^/s0 o + 0C0 Q ) + (x 1 - 110)Q o P + (z x - 40)Q Q r 


and 
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a = standard deviation of wind 
v = mean value of wind 

x = Xj coordinate of vehicle center of mass 

~z^ = coordinate of vehicle center of mass 

The response vector is related to the states, control and mean wind via the 
following equation: 

r = HjX + DjU + D 2 v (C 2 ) 

The structure of the state equations and response equations are displayed in 
Table Cl and Table C2, respectively. The values of elements of the matrices 
A, B]_, B 2 , B 3 , Hj, Dj, D 2 are exhibited from 0 to 210 seconds at five-second 
intervals in Tables C3 through C 8 . Elements which are equal to 0 will not be 
displayed, and the elements which are constant are displayed separately first. 
Numerical data was provided in English units. To reduce errors these units 
were maintained. 


a 4,l = 

1.0 

a 5, 2 

1.0 

a 6 , 3 

1.0 

II 

D— 

* 

CO 

cti 

1.0 

a 10, 9 

= 1.0 

a 12, 11 

= 1.0 

a 14, 13 

= 1.0 

a 15, 15 

= '31.6 (1/ sec) 

a 16, 16 

= -31 . 6 (1/sec) 

a 17, 17 

= - 10.0 ( 1 /sec) 
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'15,1 

1 

*16,2 

1 

>17,3 

1 

>1,4 - 
1 

a 2,l 

1 

1 3, 6 
1 

1 4, 3 
1 

1 5, 15 

1 6, 15 
J 

n 7, 16 
h l 

n 8, 16 
h 1 

n 9, 17 
.1 

n 10, 17 


31.6 (1/sec) 
31.6 (1/sec) 
10.0 (1/sec) 


1.0 

1.0 

1.0 

1.0 

1.0 

-31. 6 (1/sec) 


1.0 

-31. 6 (1/sec) 


1.0 


-10.0 (1/sec) 



31.6 (1/sec) 


d 8,2 = 31.6 (1/sec) 
djo 3 = 10. 0 (1/ sec) 
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Table C3. A Matrix (continued) 
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Table C3. A Matrix (continued) 
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Table C 3. A Matrix (continued) 
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Table C3. A Matrix (continued) 
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Table C3. A Matrix (continued) 
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Table C 3. A Matrix (continued) 
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Table C3. A Matrix (continued) 
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Table C3. A Matrix (continued) 
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Table C3. A Matrix (continued) 
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Table C3. A Matrix (continued) 
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Table C3. A Matrix (continued) 
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Table C3. A Matrix (continued) 
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Table C 4. B2 Matrix 
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0000000000000000000 00000000 < 
UUUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUIUJUUJUJUJUJUJUJUJUJUJUJUJUJUJWWUUJUJUJUJUI 

srvnsDsDsCvD'£)sCvDvCsDNDvDs 0 s£:^vC'£(% 0 s 0 '£i'£>^C'£'£iv 0 vCvC'C'£O'£)'£>v£)^'Cvr'C'C'£>'C'CvC 

— r-r—r-f— r-r-r--i — p-r-r--r— 1 — r-r— r-i — 1 — 

nnnnnnnnonnnnnrmnnnnnnnfonnnnnnnnnnnnnnnnnnn 


1 1 1 1 1 1 1 1 1 t 1 


111 1 1 1 1 1 


1 1 1 1 1 1 1 


ojrucvrvjnjruojojojrvjojrvjojouojruojruojrunjojpjrurvjnjnjojojojojojojojnjojojruojojojrvjrvj 

0000000000000000000000000000000000000000000 
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n n (*) 

O' O' & 

r- o- o- 


n n o n 
O' O' O' O' 
1 — r— r-~ r— 


n n r i i n n 

(^i Qs (J s ‘ (Js QV 
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Table C5. B3 Matrix 


TIME 

B3(?l» 1) 

B3(22* 1) 

B3(24, 2) 


1/SEC 

1/SEC 

1/SEC 

0. 

0. 

0. 

0. 

5. 

-.6879E-05 

.9844E-01 

• 1 886E-0 1 

10. 

-.98Q5E-05 

.1416E+00 

.2759E-01 

15. 

-. 12?9E-04 

• 1 759E + 00 

.3442E-01 

20. 

-.1419E-04 

.2060E+00 

.4227E-01 

25. 

-. 1633E-04 

.2337E*00 

.5095E-01 

30. 

-.1815E-04 

.2596E+00 

. 6033E-0 1 

35. 

-.19R6E-04 

. ?84lE*00 

.7209E-01 

40. 

-.2149E-04 

• 3075E* 00 

• 8485E-0 1 

45. 

-.2304E-04 

.3296E+00 

. 1028E+00 

50. 

-.24S1E-04 

.3508E+00 

. 1 225E+00 

55. 

-.25P5E-04 

,3713E*00 

•1512E*00 

60. 

-.2731E-04 

.3908E*00 

• 1861E+00 

65. 

-.28R0E-04 

,4077E*00 

,2302E*00 

70. 

-.29R3E-04 

.4225E*00 

•2858E+00 

75. 

-.3043E-04 

.4355E*00 

.3320E>00 

BO. 

-.31?4E-04 

.4470E+00 

,3399E>00 

85. 

- . 3 1 Q7E-04 

,4575E*00 

• 3207E+00 

90. 

-.3262E-04 

.4668E+00 

.2733E+00 

95. 

-.3322E-04 

.4754E+00 

•2572E+00 

100. 

-.3376E-04 

.4831E+00 

•2461E+00 

105. 

-.34P3E-04 

.4898E+00 

•2359E+00 

110. 

-.3461E-04 

.4953E*00 

.2245E>00 

115. 

-.34R9E-04 

.4993E*00 

•2158E+00 

120. 

-.3508E-04 

,50?0E*00 

. 2054E ♦00 

125. 

-.3522E-04 

.5040E*00 

. 1949E+00 

130. 

-.3530E-04 

.5050E*00 

.1914E+00 

135. 

-.3531E-04 

.5053E+00 

. 1871E + 00 

140. 

- . 3528E-04 

.5048E+00 

. 1820E+00 

145. 

- . 35?0E"04 

.5037E*00 

• 1 759E + 00 

150. 

-.3507E-04 

.5019E+00 

• 1687E + 00 

155. 

-.34Q2E-04 

.4996E+00 

• 1605E+00 

160. 

- • 3474E-04 

.4971E+00 

.1518E*00 

165. 

-.34R2E-04 

.4939E+00 

• 1419E + 00 

170. 

- . 34?4E“04 

.4900E+00 

. 1312E+00 

175. 

-.3392E-04 

.4853E+00 

. 1 196E*00 

180. 

- . 33S6E-04 

.4802E*00 

.1072E+00 

185. 

-.3317E-04 

• 4746E ♦ 0 0 

.9438E-01 

190. 

-.3276E-04 

.4687E+00 

.8091E-01 

195. 

-.3233E-04 

.4626E+00 

.6718E-01 

200. 

-.3190E-04 

,4564E*00 

.5289E-01 

205. 

-.3146E-04 

.4502E+00 

.3843E-01 

210. 

-.3104E-04 

• 4441 E ♦ 0 0 

.2355E-01 
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Table C 6. D1 Matrix 


TIME 

D 1 ( 1 4 ♦ 1) 


01(14* 2 ) 

Ul ( 14 , J) 


ft/sec 3 -rad 


ft/sec 3 -rad 

ft/sec 3 -rad 

0 . 

.2224E+04 

- 

. 17166+05 

0. 

D » 

. 2231 E + OA 

- 

. 1 7 / BE + 05 

. 2 7 J 7 F + 0 1 

1 0 . 

.2224E*04 

- 

. 1 8?6t +05 

. llo?E*02 

15* 

.2219E+04 

- 

• 1 9 a 4 £ + 0 5 

. 2 7 0 0 E + 0 2 

20. 

.2224E+04 

- 

. 19*66+05 

.4962E+02 

25. 

.2224E+04 

- 

• 21a 16 + 05 

. 7^39t+02 

30. 

.2228E+04 

- 

• 2 2 -> 2 £ + 0 5 

. 1 1 54E + 03 

35. 

.2243E+04 

- 

. 24s9E '►ub 

. 1592E*U3 

40 . 

• 2246E + 04 

- 

.26. >96 + 05 

.20876*03 

4 d • 

. 2255E + 04 

- 

.281 86 + 05 

o 2 5 6 2 6 ♦ 0 3 

50. 

.2253E+04 

- 

. 30066*05 

.30446+03 

55 • 

.2277E+04 

- 

• 3 1 a 7 6 + 0 5 

. 35896 + 0 3 

60. 

. 2284E ♦ 04 

- 

. 3 3 5 7 E * 0 5 

. 39406 + 03 

6 5. 

.2306E+04 

- 

• 35uOE + 05 

. 34396*03 

70. 

.2314E+04 

- 

.36 ) OE + Ob 

.41816*03 

75. 

. 23 1 2E* 04 

- 

.37056*05 

• 4 4 7 b E * 0 3 

80. 

.2310E+04 

- 

• 38096 *05 

• 4 7 0 b E * 0 3 

85. 

.2336E+04 

- 

.39 / 1 6 * Ob 

•63546+03 

90. 

.2332E+04 

- 

• 4 1 s2£ + 05 

. 62456*03 

95. 

.2 325E*04 

- 

•43a3c+05 

.42376*0 3 

100. 

. 23 1 9E + 04 

- 

• 46 1 26 * 05 

. 24bbE*0 3 

10 5. 

.2299E+04 

- 

• 48 0 7t + 0 b 

• 16/06*03 

110. 

.2290E+04 

- 

• 4 9 I 36*05 

.14156*03 

115. 

.2269E+04 

- 

.51 !. 4 1 * U b 

• 12616*0 3 

120 . 

. 225 1 E + 04 

- 

•52846*05 

. 96126*02 

125. 

•2225E+04 

- 

.53S9l*05 

•61566*02 

130. 

.2213E+04 

- 

.54 U6*05 

• 3206t*02 

1 3 5 • 

.2180E+04 

- 

. 54 a 96 * 05 

• 15676+02 

140. 

.2151E+04 

- 

• 54 ?s 26*05 

.9421 F +01 

145. 

• 2 1 1 8E 04 

- 

.55086*05 

. 78086*01 

150. 

.2086E+04 

- 

.55566*05 

.71316+01 

155. 

.20S3E*04 

- 

.56] 16*05 

•64106+01 

160. 

. 20 1 5E* 04 


•5658E+05 

.55106*01 

165. 

. 1897E+04 

- 

• 54^36 * 05 

• 442 1 6 * 0 i 

170. 

. 1754E*04 

- 

• 52 J 8L* 05 

• 3Sb6E + 0 1 

175. 

. 1635E+04 

- 

.50 U6 + 05 

.30546+01 

180. 

. 1540E+04 

- 

.48566*05 

• 25406 ♦ 0 1 

185. 

. 1394E + 04 

- 

• 4b 026 * 05 

.21356+01 

190. 

. 1300E + 04 

- 

.45 166*05 

.16426+01 

195. 

. 1 174E + 04 

• 

.43006*05 

. 10326 + 01 

200. 

.1077E+04 

• 

• 4 1 a 9 E * 0 5 

.13166+01 

20 5. 

•9522E + 03 

- 

• 39576*05 

• 1 1 966 + 01 

210. 

• 8587E ♦ 03 

- 

• 38 396 * Ub 

.904/6+00 
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Table C7. D2 Matrix 


inE 

02(11* 1) 

02 (J 2* 1) 

02 l l 3* 1 ) 


LB-SEC-RAD/FT 3 

LB-RAD/FT 3 

1/SEC 

0 . 

0 . 

0 • 

- . 5 0 5 1 E - 0 2 

5. 

• 5 8 4 0 L “ 0 i 

. 1 0076-03 

- . 4 0 n / L - 0 2 

10. 

. 1 1956+00 

. 2 0 b 7 K " 0 3 

-.3 Is 16-0 2 

15. 

•lbjOE+0 0 

.31946-03 

- • 2 <2 0 7 E - 0 2 

20. 

• 2 ^ 2 3 E ♦ u 0 

• 4 2 5 b E — 0 3 

- . [2996-02 

25 • 

•30 i2E+uO 

•52lbE-UJ 

- • 3 7 9 9 E - 0 J 

30. 

.35 Wti + uO 

. 6 U b 8 E - 0 3 

• 4 3 4 0 £ - 0 J 

35. 

•4Q62E+00 

. 6 7 3 56-0 3 

. 1 0 9 d E - 0 2 

40 . 

• 44 /86+uO 

. 7 1 2 5 E “ 0 J 

. 1 6 '4 .2 E - 0 2 

45. 

.48^46+uO 

. 7 3d 7 £-03 

• 2 0b 7 E-0 2 

50 . 

. 5 0 4 6 1 . + 0 0 

. 7 9 4 2 E - 0 3 

• 2 4 8 2 E — 0 2 

55. 

. S2506 + U G 

. 7922E-G 3 

• 2 8 i 6 E - 0 2 

60 . 

.53086+U0 

. 7707E-03 

. 2 6 6 3 £ - 0 2 

55 • 

• 5 2 4 0 E + 0 0 

. 7 4 « & t - 0 3 

. 24Y3E-U8 

70. 

.50 SOti + OG 

. 7 70 It- 0 3 

. ? b ^ 5 £ - 0 2 

7 5. 

• 47m 7t + 00 

.44366-03 

. [ bsu£-G2 

80 . 

.43 5 5 £ + O 0 

. 31 3 2 L - 0 3 

• 1 3 3 6 L - 0 2 

8 5. 

•38746+00 

. 20 94E-0 3 

. { 0 4 7 E - 0 2 

90. 

• 3 3 1 66 + 0 0 

. 1 30 7 E-0 3 

. 9 7 4 b £ — U J 

95. 

.2 /?96+00 

. 74bOL-U4 

. 7 1 l 0E-0b 

100. 

.21 ‘/8t + u0 

.420 6E — 0 4- 

.57 JbE-0 J 

105. 

. 1 7 <96+00 

.23796-04 

. i+ 4 4 0 £ - 0 3 

110. 

. J 3- 6t. >00 

. J 335E-04 

• 3 J 9 5 1 - 0 3 

lib. 

• 1 1 o 5 E + 0 0 

. 7 4 1 2 E - 0 5 

• ? 6 0 9 E - 0 J 

120. 

.88216-01 

.4271 E-Ob 

• 2 0 1 3 £ - 0 3 

125. 

.70156-01 

« 23996-05 

. 1 4^ JE-0 J 

130 . 

.55 ?4E-01 

.13 7 BE- 06 

. ] 0 866-0 3 

13 5. 

.4433E-01 

. 7 b 2 3 £ - 0 6 

. 7 5 6 9 E - 0 4 

140 . 

• 35 3 7ci~0 1 

„ 44b0£ — 06 

• 5 3 J 1 £ - 0 4 

145. 

.28336-01 

. 2b2 7E-0b 

• 38 <dE-04 

150 . 

.22-lL-ul 

• i 4 7 0 £ - 0 b 

. 2509£-04 

155. 

• 1 d'4.9£-u 1 

• 8 b 9 b E - 0 7 

• 1 75b£-04 

160. 

. 15 i 06-01 

.51d4£-u7 

. 1 20 7 E - 0 4 

165. 

.12926-01 

•2956E-0 f 

• 7 2 b £ - U b 

170. 

. 1 0296-0 1 

. 1 772E-0 7 

• 4 7 5 b E - 0 5 

175. 

• 86 1 46-02 

. 10 72E-07 

. P904E-05 

180. 

• 73*)b6 — 02 

. b 2 7 7 E — 0 6 

.11 o6£-03 

185. 

.62o9t-02 

.3339 E - 0 8 

- • 299 3 £-06 

190 . 

• 52 9 7t.-02 

• 1 5 / 1 E-0 8 

- . 1 3 -49£-0b 

195. 

•45o6E~o2 

.4U58E-09 

- • 2 2 7 E - 0 5 

200 . 

• 38456- o 2 

- . 30856-09 

- . 3 0 b 8 E - 0 5 

20 5 * 

.32-36-02 

7217E-09 

- . 3b £-05 

210. 

• 28046-02 

- . 8 7 J7E-09 

-.39 9 / E - 0 d 


L)£ ( i A * 1) 

1/SEC 2 

-.44b2E~02 
. 198/E-U1 
. 4 2 5 / E - 0 1 

• b40b£-0 1 
• 79E-0 i 

. i u n t ♦ o o 

.11 / J E + 0 0 
. 1 3 2 b E ♦ 0 0 
. 1 4 .3 1 E> U 0 
. Ib'-oE + OO 
. 1 o 9 9 E + 0 0 
• iolOE+00 
. 1 d 9 4 E ♦ 0 0 
. I9/3E+00 
. 2054E+ 00 
. i b 3 2 E ♦ 0 0 
. 12/2E+00 

• 1 U 2 0 E + 0 0 
. 7ti4bE-Ul 
.571 7E-01 
.41 1 JE-01 
. 29 J9E.-0 1 
. 2u94E~0 l 
. 1 /24t-0 1 

• !39bE-Gl 
. 1 lObE-Oi 
. 798 IE-02 

.o/9/£-02 
•3900E-08 
.233bE-02 
. i b 69£- 0 2 
.5 774E-03 
.421 / E - 0 3 
•42G/E-03 

• 2 2 4 2 E - 0 3 
. i ( 42E-U 3 
. 3279E-03 

• 5436E-0 3 
. ^0 /E-0 3 
. 9 b 7 b E - 0 3 
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Table C8. HI Matrix (continued) 
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Table C8. HI Matrix (continued) 
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Table C8. HI Matrix (continued) 


i — ! "vj .\j 'vj ?\j r\j -vj\i <\J '\j "vj ""i "'0 -o nn m ^ ^22 2 

JOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO® 

j uJ Ixl UJ UJ UJ uJ UJ 


■ O T> 

■ o O' : 

i -j "if 


'vy.Mnnn-iOBM.ftH3 
, , „ _ . . j o iM ^ B s o n-n n o m ^ n 

, < rO '''"14 4 O' 4 


' n n -< o ,m t xi ® o n 

V CnBBM-'BBMBT' -.M-iO ft- | BB7'M' 

• n n h n n o n o /» o o 
. ,-*n _ ry <n j- .n h ^ n n o' - 


« -h tv <v: 'vj !\, f\j oj ^ n m n^nnnnnnwM rvi 

OOOO =>30 0 0 000 oooooooooo o 

♦ ♦♦ + > + + + + + + + + + + + + + + + ' f + + 

lx! UJ UJ UJ UJ UJ Ul U.I UJ UJ uJ Ixl UJ LU Ul i.U U.J LiJ UJ UJ UJ UJ UJ 

4- o o 3 h m j* ^ H c 3 o.^.-.-«'n.T)-*o.n'|- 

^ N. } O O ft ^ ^ <ti O O »*) '*M 3 -S’ ift M 3 4 O ft) 

r nT>to n ■" — » ~) ft).^ 00034 o 3 4 m x> 3 

n h m ,t n s t> - • — • — • — • O' >“ ^ 


>0 000 
J UJ UJ UJ llJ 


ooooooooooo 


O ftVI <*> ■ 

o n "o 

4 V — < 


UJ IjJ UJ 
t> 7* O 
— ■ 3 J; 


UJ UJ UJ UJ UJ 

o o n m n 

f — 4 ft 4 JO 

-o 'O 4- -• > 

TJ -O 3 CVi 4 r 


Jj LiJ UJ 
— < 4 'O 
IfKft - 
o 'O 


- o 

4 UJ 

to 


2 Saooooo 30 = 3 = SSS 


o o o o o o 


> 000000000 ; 


UJ UJ ; xl ixl UJ ! xl uJ Ul UJ 'xl *jj i-d UJ UJ -iJ UJ uJ uJ UJ >3 -Jj UJ UJ :*J UJ UJ xJ UJ uu 5 
n ji j <> o J’ M v £■ O O n X 4 -T> J”> >* — • AJ 3 3 IT) -f- 4 ^ •£ 0^ • ^ ^ 

n!ooJ oyfi - or oVMoo^Mnt^'OE^B'Jft- 

fM (V n n j J If if: ft- ft ■ ^ s , — . — j — r— r— — r- r~ -.r a) a. x x. <x 


; >jj >jj UJ UJ -xl ‘xl UJ ijJ UJ uJ UJ 
■ n j \i -i o m J M .ft 


;y i\j c\j M m n (*) <n r> n n n n n n n ro m n n n 4) cvj rvj *v -v -\j oj <\j aj v aj n «v ^ ^ 2J JJ JJ JJ J; 

oooooooooooooooooooooooooooooooooot^oooooooo 

i tL tL UJ UJ U UJ UJ UJ it) UJ UJ UJ i UJ lL UJ UJ UJ Ul UJ ui UJ UJ i uJ UJ UJ UJ Ixl UJ LlJ 13 uj Ul UJ bJ UJ bJ uJ UJ UJ U 

5 SJS-SS25SS22S2S2S2J?SSSS?5nSi?S-XSS22£c&S^32| 


f'linninnntnnnnnnnnrtrtinnnnnnn.nnnnnn 
OOOO o O O O O O O O O O O O 0300303000000 


♦ + 


J UJ UJ 
■ r» ro 
- a; 


til UJ UJ uJ UJ UJ UJ UJ IjJ i4J 'xJ IjJ UJ ’jJ ’jJ ‘jJ UJ UJ UJ !jJ UJ 'jJ UJ tiJ UJ UJ UJ UJ n 

c\’ •£ o <1 vT m r- o 'ft a o «— * m ’J* ift ir ^ ^ o ^ ■ o vft ^ 'C?' 1 

a: w n n n j- ftftijifti’ftftiftiftiftftft'Oftftiftiftft'ft'ft ^ n n n m im ' 


: 'jJ txl '-J 
IT ^ O' ' 
a ft ft ^ 1 


i i r i i l l l i l ■ > 


i i l i I * 


-h iiM m <vj (v t*> r> *> •?*> m 'no -no no r) m »*> m no r> n m n no m n n ro no n n n •<*> n m n fto 

0000000000000000003000000000000000000000000 


♦ ■+ '♦ . 


«♦ • ♦ ♦ 


•♦♦♦♦♦♦•♦ 


+ ♦ 


♦ ♦ 


♦ + ♦ ♦ 


+ '♦ 


• ♦ ♦ 


-♦ + 


UJUJUJUJUJUJUJUJUJUJUJUJUJUJuJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJuJUJUJUJUJUJUJUJUJUJUJUJUJ 

ao ft> - t®ftOMU'ft(?'lftOM\)ftOOOONOtJ'noMft3'N- 

.No-<nNi/i4'N(MOco(j'4(\Jftffl4-ioO'-<niftcDonN<-<ift 

V £ )r r ) r*. l ^rvjo4'oor r )Otfi^r~mr^— innHC04'ON'tftjoa'®Nftft>ftft'0'0ftft'0^f^ ,v< D® 

nHftjfts-i-iHMMn'tftftitftftftftftumftift^ftj'ftnnnnnnnnfftnnrftnnnnn 


) *n r> r»> -*v-»vo rn -n ■»> *o -m "n ^ *o r> ro ■■ 

SOOOOOOOOOOOOOOOOOOOO( 


♦ ♦ ♦ <♦ 


UJ UJ UJ UJ 
O' <t LD \0 n 

niftftftn 
o<rrt\jr- 
— • rg .^n lt> ft) 


UJ UJ UJ 
so r- <r 
oj oo 4- 


UJ UJ UJ UJ L 
ft n oin ( 
^00*^1 
n u-> r» o - 


♦ ♦ ♦ ■ 


♦ ♦ ♦ ■ ♦ 


,^^^-w^.UJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJ 
ift)j)ftHft)(\io-<Nrtfttcft)ocoMftO'nNsO'a'®'-<no-i 
)XO4fflMft)Mri-H4n0'4C0 00 X3N4'-'ft)o4' N (MMinNa' 

jftjnnn>t4'tift , /uft>Oft'ft'Oft-M^ f * ,N ®® , H®®®®®^ 


I 


< i 


I I I I I I I I I 


I i i - 1 i i I I 


i i i i 


i i i 


i i i i i 


UJ UJ i 
-I sO 3 
t/> 'NJ t 


■ rvi ift> ® — < -h . 


J UJ UJ UJ 
» ft) in r\J 

) -4 ^ sO 

■ r- O' — < 


O' 0J 4 - 
on jr 
n <f in l 
oj ru rvj r 


♦ ♦ ♦ 

J UJ UJ UJ 
) — « O' ao 

- lT> r- 

lOMft 

j rvj 


» ill 

UJ 


UJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJUJ 
moNBft-<-'ftftftJNft)ftH\J[ r >ftnojiBo 
N,n^®NN-'®^t)COM-<ION-Hl/)JlOO 

fMo-.f\jr-i^O'-.sD^inrvjin^4-f-nao^ooir 

— <-hX\D 4 nrvjrvj-'^'B'O'J'nM-t-'O^fti'J 


UJ UJ 
o r\j 
o m 


0 o 

1 I 

■xl UJ 


I I I I I I I 


I I I I I I I I 


I I I I I I I I I I I I I 


w to 
— < I— 

X u. 


I I ♦ + + 

IjJ UJ i4j 'jj ’xl 

n ^ ft ^ n 

ftj rvj •& o — • 

ft? 4- o 4 1 r 

rvj ^ rvi 4- n- 


> O O O 3 O C 

♦ + ♦ ■ 

IjJ UJ >xl l 

o ft) ft) r 

n m in r 

— « ru m : . 


J UJ UJ 'aJ 
) o ft) ID 

• ON® 

i ft) -• m 


UJ UJ UJ '. 

O O O 

CO 4 1 N 

ft) X ft) V 


. ox -i -h -i -h . 


1 UJ UJ 'xJ UJ UJ UJ UJ UJ f 4j : xl UJ UJ «xJ UJ UJ Ixl UJ ’xl UJ UJ UJ UJ 'xl 

i 4 jft'ftjft J’T'ftftftMB-'nftnBonoo-tft. 

,-tr-o4!BON4nft54ft)0'(\in(\J4’7'BincOftJN: 
ift)r*->!\jOf'~ft>U‘)44 , ft)X) r Ur'~4‘'- <— •rOO'7'^H4 3'4< 
,^,-H^co^ft)ft)4ftorvjf\j-H^— •o'r^'ft)4-4'nf\i(Mft 


I I I 


I I 


I I I I I 


I I I I I I I I I 


I 


I ♦ ♦ ♦ 

uj uj uj uj uj 

N O ft) D 
4 ft) 3 ft) 3 
3 C\J 4 O O 
j h m ft) n 


>0000000330000'; 


i OOJO ^ M M n M IM M M M M M M M 

>00000003 0 000000000000 

I I I l l I l l l t I l l i I I i 


UJ Ixj 

ru o 

O 4 


UJ UJ UJ 
■ft) O' 4 
o 3 m 
4 BM 
-« -h rvj 


UJ UJ UJ 
O X> XT) 
3^0 

n n o' 
rvj oj nj 


UJ UJ UJ L 
o O' O' t 
-*44; 
3 3 ~| r 
rvj rvi m - 


jijJUJuJUJUJt3UJUJUJUJUJUJUJuJUJl3UJUJi3UJUJUJUJUJ 
)ftft0'0'ft)4--iMMOft)ft)'inftft'-' H 0'4Hft(\J!ft 
< n n n ft) n O' o ft) o "" ft) m ft ^ m — * \i .T' j' 3 

’ 04-^4 ftcoft) 33 40-H.^)ft~03'30— »;\J— I J' N- 3 

<-^x3nr\J— <3— *3— 'rvjrommnrvjrvjoj 


i i 


i i i i 


i i 




73 


210 



Table C8. HI Matrix (continued) 
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Table C8. HI Matrix (concluded) 


TIME 

HI ( 14.24) 
FT/SEC 3 

0. 

0 . 

5 • 

. 1 bbbE + 01 

1 0 . 

• 688 7 E ♦ 0 1 

15. 

.lb7u£+02 

20. 

.2 784t+02 

25. 

• 42bbE+02 

30. 

.bob/L+0 2 

35. 

. 7<+62£ + 02 

4-0 . 

. 89 1 HE + 0 2 

4-5 • 

. 1008E+03 

50. 

. 1 1 09E+03 

55. 

. i iblE+03 

60. 

. 1 1 8 1 E ♦ 0 3 

6 5. 

. 1209E+03 

70. 

. 1 336E + 0 3 

75. 

. 1040L+03 

80. 

.91 7 i t + 0 2 

85. 

. 722b£+02 

90. 

.42bvE+02 

95. 

. 55bbE + 01 

100. 

-.2E>2bE + 02 

105. 

- ,4b43ti + 02 

110. 

- . 639 Jt + 02 

115. 

- . HhOhE+02 

120 . 

-,932bE+02 

125. 

- . 960 8E + 02 

130. 

- .9398E+02 

135. 

-.9040E+02 

140 . 

- . 8826E ♦ 02 

145. 

7944E+02 

150. 

74lb£+02 

155. 

-.687 bt + 02 

160 . 

- • 6368E + 02 

165. 

- . 590 7E + 02 

170. 

- ,b442£ + 02 

175. 

-.602 7 £+02 

180. 

- . 468bE + 02 

185. 

- . 4339E ♦ 02 

190. 

- . 40 l+fE + 02 

195. 

-.3700E+02 

200 . 

-.341 4t + 02 

205. 

- . 3 1 4GE + 02 

210. 

-.2890E+02 
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APPENDIX D 
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APPENDIX D 
TRAJECTORY DATA 


The nominal trajectory is described quantitatively in Table Dl. 


77 



Table Dl. Trajectory Data 
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Table Dl. Trajectory Data (concluded) 
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